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Abstract 

Using full-potential density functional theory (DFT) calculations, we found a small asymmetry 
in the Fermi velocity of electrons and holes in graphene. These Fermi velocity values and their 
average were found to decrease with increasing in-plane homogeneous biaxial strain; the variation 
in Fermi velocity is quadratic in strain. The results, which can be verified by Landau level 
spectroscopy and quantum capacitance measurements of bi-axially strained graphene, promise 
potential applications in graphene based straintronics and flexible electronics. 
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1 Introduction 

Graphene, a single layer of carbon atoms ar- 
ranged in a two-dimensional (2D) hexagonal lat- 
tice, displays remarkable mechanical and elec- 
tronic properties promising many applications 
in nano-devices [jj. The novel properties of 
graphene arise from the linear energy dispersion 
near the K point of the hexagonal Brillouin zone 
(BZ): 

E-EF^±VFhk (1) 

where Ep is the Fermi energy, vp is the Fermi ve- 
locity of electrons/holes near the K point of the 
BZ and {hk) is the momentum. Many graphene 
device characteristics depend on vp. For in- 
stance, graphene's fine structure constant 



(S. I. units) 



appears in the graphene device characteristics of 
graphene field effect transistor (GFET) [2j; for 
an ideal graphene with a uniform channel poten- 
tial Vch, the quantum capacitance of graphene is 
given by [3]: 
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where T is the temperature and other symbols 
have their usual meanings. Thus, any change 
m. vp would affect the graphene device charac- 
teristics. Although strain-induced variation in 
vp has recently been observed experimentally 
[1] by Raman spectroscopy study of uni-axially 
strained graphene, the quantitative dependence 
oivp on strain remains unclear. Here, we report 
our theoretical investigation on the effect of biax- 
ial strain on the Fermi velocity of charge carriers 
in graphene, which is experimentally known to 
sustain in-plane tensile elastic strain in excess of 



20% [5]. Thus, our study mimics the experimen- 
tal condition where graphene is supported on an 
ideal flat stretchable substrate. 

2 Computational Methods 

We use the DFT based full-potential (linearized) 
augmented plane-wave plus local orbital (FP- 
(L)APW+lo) method ^ as implemented in the 
elk code [7j- For the exchange-correlation term, 
we use the Perdew-Zunger variant of local den- 
sity approximation (LDA) [8j, the accuracy of 
which has been successfully tested in our previ- 
ous works [9]. For plane wave expansion in the 
interstitial region, we have used |G + k|maa; x 
Rmt = 9, where Rmt is the muffintin radius, 
for deciding the plane wave cut-off. The Ap- 
point grid size of 30 x 30 x 1 was used for all 
calculations. The total energy was converged 
within 2^ eV/atom. The 2D hexagonal struc- 
ture of graphene was simulated by 3D hexagonal 
super cell construction with a large value of c- 
parameter (|c| = 40 a.u.). The application of 
homogeneous in-plane biaxial 5 strain was sim- 
ulated by varying the in-plane lattice parameter 
a(= |a| = |b|);(^ = (a — ao)/ao, where oq is the 
ground state in-plane lattice constant. 

3 Results and Discussions 

In Figure 1(a), our calculated results show that 
under compressive strain of 5% graphene does 
not show any buckling, i.e. its planar structure 
remains preserved. Figure 1(b) depicts the linear 
energy band dispersion of unstrained graphene 
near the K point of the BZ. In Table 1, we com- 
pare our calculated values of and vp with 
some reported values. Our calculated variation 
of Fermi velocity of electrons vp{e), holes vp{h) 



1 



-3 




0.00 0.05 0.10 0.15 



Buckling parameter: A (A) 



(a) 




K 

(b) 



Figure 1: (a) Buckling probe of graphene at 5% 
compressive biaxial strain using the energy min- 
imization procedure. Inset shows the side and 
the top-down views of buckled graphene in which 
alternate atoms in a hexagon reside on two dif- 
ferent parallel planes; buckling parameter A is 
the perpendicular distance between these two 
planes; for planar graphene A = O.OOA. (b) En- 
ergy bands of unstrained graphene near K point 
of the BZ. 



Table 1: Calculated values of oq and vp com- 
pared with reported values. 
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Experiment, single 
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and their average value vp = [vpi^) + vpih)] /2 
at different strain values are shown in Figure 2. 
Since vp{h) values are slightly greater than vp{e) 
values, the band structure of electrons and holes 
are not exactly symmetric in qualitative agree- 
ment with a recent experiment |12j and the ori- 
gin of this asymmetry is not well understood |12j . 
The decrease in Fermi velocity with increasing 
strain is due to a reduction in the vr-orbital over- 
lap [1]. Our calculated data in Figure 2 best fit 
with the following equations: 

vp{e) = 0.82386- 1.26826 X 5+0.71355 xJ^ (4) 

vp{h) = 0.84161 - 1.3269 x (5 + 0.76501 x (5^ (5) 
VP = 0.83273 - 1.29758 x 5 + 0.73928 x 5^ (6) 

A possible precise measurement of the variation 
oi vp with 5 may be carried out with Landau 
level spectroscopy study of a bi-axially strained 
graphene, since the Landau level spectrum is 
given by [12]: 

En = Ep) + sg Ti{-a)vp^j2ehB\n\ (7) 
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e 2: Fermi velocity variation with strain 6. 



where n is an integer, is the energy at the 
Dirac point, B is the magnetic field perpen- 
dicular to graphenes plane and other symbols 
have their usual meanings. Measurement of the 
quantum capacitance [3J of bi-axially strained 
graphene can also be useful. 

4 Conclusions 

There exists a small asymmetry in the Fermi ve- 
locity of electrons and holes in graphene whose 
reason is not yet clear. The Fermi velocity in 
graphene can be modulated by strain engineer- 
ing for potential applications in graphene based 
flexible electronics. 
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